Uncertainty relations in the realm of classical dynamics 
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It is generally believed that classical regime emerges as a limiting case of quantum theory. Explor- 
ing such quantum-classical correspondences in a more transparent manner is central to the deeper 
understanding of foundational aspects and has attracted a great deal of attention - starting from 
the early days of quantum theory. While it is often highlighted that quantum to classical transition 
occurs in the limit ft — Y 0, several objections have been raised about its suitability in some physical 
contexts. Ehrenfest's theorem is another widely discussed classical limit - however, its inadequacy 
has also been pointed out in specific examples. It has been proposed that since a quantum mechan- 
ical wave function inherits an intrinsic statistical behavior, its classical limit must correspond to a 
classical ensemble - not an individual particle. This opens up the question "how would uncertainty 
relations of canonical observables compare themselves in quantum and classical realms?" In this pa- 
per we explore parallels between uncertainty relations in stationary states of quantum systems and 
that in the corresponding classical ensemble. We confine ourselves to one dimensional conservative 
systems and show, with the help of suitably defined dimensionless physical quantities, that first 
and second moments of the canonical observables match with each other in classical and quantum 
descriptions - resulting in an identical structure for uncertainty relations in both the realms. 

PACS numbers: 03.65.Sq, 03.67.Ta, 05.20.Gg 



I. INTRODUCTION 

It is imperative to retrieve classical dynamics as a lim- 
iting case - in its domain of validity - from quantum 
theory. The generally prevailing notion is that classical 
mechanics emerges in the limit h — > 0. Yet another clas- 
sical limit discussed widely is the one following Ehren- 
fest's theorem. Pedagogic discussions, in several text 
books on quantum mechanics, are essentially confined to 
these two approaches towards classical regime. However, 
both these routes are shown to be not universally satis- 
factory [ll-Q • It has been pointed out that classical realm 
- resulting from a quantum mechanical state - is ought 
to correspond to an ensemble - not a single particle [1,11] . 
The averages, variances and other higher order moments 
of the quantum and classical probability distributions are 
therefore expected to agree in the limiting case. 

In order to compare the statistical form of classical dy- 
namics with the corresponding one in quantum dynam- 
ics, phase space probability distribution of the classical 
ensemble (a counterpart of the corresponding quantum 
state) needs to be identified. The classical phase space 
probability distribution satisfies the Liouville equation 
and the phase space averages of the classical observables 
are shown to exhibit analogous dynamical behaviour as 
that of the corresponding quantum case - even when 
Ehrenfest's theorem breaks down 

Another approach, when one confines to stationary 
state solutions of the quantum Hamiltonian, is to graphi- 



cally compare the probability density function Pq^[(x) = 
\ip n (x)\ 2 with the corresponding classical probability dis- 
tribution Pcl (%) of an ensemble and to recognize that the 
envelope of the quantum probability density approaches 
the classical one in the large n limit [9] . 

In this paper, we show that the first and second mo- 
ments of suitably defined dimensionless canonical vari- 
ables evaluated in the stationary states of one dimen- 
sional conservative quantum systems match with those 
associated with the corresponding classical ensemble. 
This, in turn leads to identical structure for uncertainty 
relations of the dimensionless position and momentum 
variables in both classical and quantum domains - bring- 
ing out the underlying unity of the two formalisms - irre- 
spective of their structually different mathematical and 
conceptual nature. 



II. CLASSICAL PROBABILITY 
DISTRIBUTIONS CORRESPONDING TO 
QUANTUM MECHANICAL STATIONARY 
STATES 

We begin here by reviewing the classical probability 
distributions Q for an ensemble of particles bound in 
one dimensional potentials V(x). The probability density 
function for position of a single particle, whose initial 
position and velocities are specified, is given by 



P 



CL 



\x) = 5[x-x(t)] 



(1) 
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where x(t) denotes the deterministic trajectory of the 
particle at any instant of time t. However, the quantum 
mechanical probability density Pq'm(x) — \^p n (x)\ 2 asso- 
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ciated with the stationary state solution ij) n (x) of the sys- 
tem, is not expected to approach - in the classical realm 
- to the single particle probability density of Eq. ([I]) 
rather, the locally averaged quantum probability density 
does approximate to a probability distribution Pq~l(x) of 
a classical ensemble of particles (of fixed energy E) in the 
large n limit 0. 

The phase space probability distribution Pql(x,p) 
of an ensemble of classical particles of fixed energy 
E and bound in a potential V(x) is proportional to 



p 

2m 



V(x) - E 



obtained by integrating over the momentum variable p: 



dpP C h(x,p) 



Constant ■ 



dp 8 



•(2) 



Using the properties 6 (ax) = ^S(x) and 5(x 2 — a 2 ) — 
w- [S(x + a) + 8(x — a)] of the Dirac delta function, the 



. The position probability function is classical probaility distribution reduces to 



P CL (x) = Constant • J dp2mS (p 2 + 2m[V(x) — E] 




2m 



V / E-V(x) 7 



[E-W)]j dp HP+V2m[E-V(x): 



+ 6 (p - y/2m[E- V(xj 



(3) 



where M denotes the normalization factor, such that 



J- A dx Pql(x) — 1 (here the integration is taken between si° n Pcl(x) — \ 



the classical ensemble reduces to the well-known expres- 



Va 1 



(—A, A) as the probability distribution Pcl(^) vanishes 
outside the classical turning points |x| > A). It may be 
readily seen that, by substituting E = ^muj 2 A 2 , and 
V(x) = ^muj 2 x 2 , in the familiar example of harmonic 
oscillator, the position probability distribution Eq. of 



The phase space averages of any arbitary function 
F(x,p) of position and momentum variables get reduced 
to those evaluated with position probability distribution 
function Pcl(x) as follows: 



(F{x,p)) C L = 



dx J dpP C h(x,p) F(x,p) 
Constant J dx J dp 8 f^— + V(x) - E^j F(x,p) 

E _ 2 y (a;) J dp [d (p + ^2m[E-V(x)]) + 6 (p - y/2m[E -V(x)])] F(x,p) 



Constant • / dx 



Constant ■ / dx 



2m 



E - V(x) 



F(x, -y/2m[E-V(x)]) + F(x, y/2m[E - V(x)]) 



dxPch(x) F (x, -y / 2m[E - V(x)] + F (x, y / 2m[E - V(x) 



(4) 



We define dimensionless (scaled) position and momen- 
tum variables as follows, 

such that \X\, \P\ < 1 in a bounded system. 



In the next section, we compute the first and second 
moments (X) C l, (X 2 ) C l, (P)cl, (P 2 )cl of the classical 
probability distribution in three specific examples of one- 
dimensional bound systems. We then compare these 
classical averages with the quantum expectation values 
(X)qm, (X 2 )q M , (P)qm, (P )qm - evaluated in the sta- 
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tionary states ip n (x) and identify that they agree with 
each other in the classical limit. 



III. COMPARISION OF FIRST AND SECOND 
MOMENTS OF THE CLASSICAL 
DISTRIBUTION WITH THE STATIONARY 
STATE QUANTUM MOMENTS 

We focus now on three specific examples of one dimen- 
sional bound systems viz., harmonic oscillator, infinite 
well and the bouncing ball - both in the classical and 
quantum domain. We evaluate first and second moments 
of dimensionless position and momentum variables (see 
Eq. for definition) and identify that the quantum mo- 
ments - evaluated in stationary eigen states of the Hamil- 
tonian - match with their classical counterparts. 



A. One dimensional harmonic oscillator 

The classical probability density for finding a system 
of harmonic oscillators - all having same amplitude A 
- between position x and x + dx is given by (see the 
paragraph following Eq. ((3])) 

CL[) \=0, for|.T|>A [) 

We consider scaled canonical variables X = j| and 
P = -, p = = — —, r and evaluate the averages of 

X, X 2 , P, P 2 (by making use of Q and ©) as given 
below : 



(X) 



CL 



(X 2 ) 



(P) 



(P 2 ) 



CL 



CL 



CL 



I /^cl(x)* = ^ f_ A dx 

i /^^^i I_ A dx 
l — J A A dxP CL (x) LJ^~ 



x 



\[A^ - X 



2 



x~ 



V 'A 2 - x 2 



I 

2' 



2 muj . 



E mu 2 x 2 

2 



'2m 



E raw 2 x 2 

2 



i 2 oj 2 A 2 



dx Pcl (x) 2m 



-A 



E mux 2 x 2 

2 



A 2 IT 



dx \/A 2 



(7) 

(8) 

(9) 
(10) 



-A 



r 



The variances of X, P are given by, 

(AX) 2 CL = (X 2 ) Ch -(X) 2 CL = ± 

(AP) 2 CL = (P 2 )cl-(P) 2 cl = \ (11) 

and hence the product of variances obey, 

(AX) 2 CL (AP) 2 CL ^± (12) 

in a classical ensemble (characterized by the probability 
distribution ©) of harmonic oscillators. 

The stationary state solutions of the quantum Hamil- 
tonian 



are given by, 



tj P 2 _l 1 2-2 
H = — mcu x 

2m 2 



1/2 



(13) 



^ )= nSr) H n {^Thx)e-^ (14) 



where H n are Hermite polynomials of degree n; the cor- 
responding energy eigen values are, 



n + — ) iu: . n 0. 1.2. 



(15) 



The classical turning points associated with the energy 
E n are readily identified as, 



A n 



2E n _ (2n + l)h 



(16) 



We consider scaled position and momentum operators, 



X = — = x 



p = 



(2n + l)h 



y/2mE n yj(2n + l)hi 



(17) 



_ p_ 



corresponding to the classical ones X — 4, P 

The expectation values of X, X 2 ,P, and P 2 , evaluated 
in the stationary states 4>n(x) (see Eq. (fT4")l ) are given by, 
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(X)qm 



(X*)qm 



(P) 



QM 



(P 2 )qm 



(2n + 1) h J_ c 
mix) 



dx \^ n (x)\ 2 x 



h2 n n\ ^{2n+l)it 
muj 



dx ( yj mix) /hx) e h x = 0, 



(2n + l)fi 

virJ 



|V> n (x)| 2 x 2 
1 



2 n n!(2n+ 1)^ 



y (2n + 1) muj 

—i 

2™ n! ^(271 + l)7r J- 

-h 



dx P 2 ( \J muj /hx) e R x 2 = 

o 

dlp n (x) 



dx 



r 2 d 



dx H n (y/mLu7hx) e h — ( H n (\/ mixs/hx) e 

dx V 



0, 



(2n + l)mw 



dx?/>*(x) 



d 2 ^„(x) 

dx 2 



1 



mw7r 2™n! (2n + 1) 



dx P n (y 'muj /Hx) e w 



dx 2 



H n (y/m^J/kx)e- Jm ^) = -. (21) 



(18) 



(19) 



(20) 



Clearly, the quantum expectation values Eas.(|18 |) -(|2i p 
match with the classical ones given in Eas.([7|)- ([TU| and 
we obtain the uncertainty product - for all stationary 
state solutions of the quantum oscillator, 



AX 



QM 



AP 



1 

QM 4' 



(22) 



by, 



B. One dimensional infinite potential box 



We consider a symmetric infinite potential well defined 



V(x) 



for - f < x < § , 
oo for |x| > -|. 



(25) 



It is pertinent to point out here that the commutator 
relation, 



[X,P] = 



P 



1 i 



(2n+l)h 



[x,pl = 



2n + l 



leads to the uncertainty relation, 



(23) 



The particles move with a constant velocity within the 
box and get reflected back and forth (i.e., the momen- 
tum changes sign when the particles hit the walls). The 
position probability distribution for an ensemble of clas- 
sical particles confined to move within the box (such that 
— ^ < x < -|) is a constant (as can be readily seen by 
substituting Eq. (|23|) in Eq. Q and is given by [9J, 



^cl(x) 



for jxj > | 



(26) 



(axY (a p] 2 > — — — -. (24) 

V Jqm V Jqm ~ 4(2n + l) 2 v ' 

In the large n limit one obtains the result 
lim^oo (AX J AP > - which is usu- 

V / QM \ / QM 

ally expected in the classical regime. However, the exact 
result (|2"2l for uncertainty product holds for all the 
stationary state solutions - and strikingly, this result 
matches with that of a classical ensemble of oscillators 
with fixed energy E (see l|12p). 



which obeys J_ L / 2 Pch(x) dx — 1. 

In this example, the dimensionless position and mo- 
mentum variables are identified as, 

X = X P = P = JL (27\ 
(L/2)' V2^P |p| 1 ' 

and the classical averages (X) C l, (X 2 )cl, (P)cl, (P 2 )cl 
are readily evaluated using the probability distribution 
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(X) 



CL 



2 

Z 2 



L/2 



L/2 



cdx = 0, 



dxP C h(x) 



L 2 /4 



(P) 



CL 



4 
0, 



L/2 



-L/2 
(^ 2 >CL 



1 

3' 



1. 



(28) 

(29) 
(30) 



So, we obtain the variances of X, P as, (AX)^, L = | 
and (AP)q L = 1 for the classical ensemble of particles of 
fixed energy E, confined within the infinite well, which 
results in the following variance product: 



(AX) 2 CL (AP) 2 CL = -. 



(31) 



Quantum mechanical stationary state solutions (even 
and odd parity) for a particle confined in an one dimen- 
sional infinite potential well (1251) are given by, 



(+)/ 



— cos(n7r x/L), n — 1, 3, 5, . 



V4 >{x) = \j — sm(mrx/L), n=2,4,6,. 



and the corresponding energy eigen values are, 

n 2 



En — 



2 m L 2 



(32) 



(33) 



The scaled dimensionless position and momentum oper- 
ators (analogous to the classical ones (|27p) are chosen 
as, 

X = — P= t = P (341 
L/2' y/2mE n tiitH/ L ' K 1 

The expectation values of X, X 2 ,P, and P 2 are eval- 
uated in the stationary states (both even and odd) to 
obtain, 



(X) 



(x 2 : 



(p, 



QM 



QM 



QM 



-i- / dx\^+/-\x)\ 2 x = 

L / Z J-L/2 

dx\^+'-\x)\ 2 X 2 



(35) 



1 


,L/2 


L 2 /4 


J-L/2 


1 


2 


3 n 


2 7T 2 ' 


. L 


,L/2 


nix 


J-L/2 


o, 





(36) 



dxi>W-\x 



d<ip { +'-\x) 



(P 2 ) 



L' 2 



QM 



L/2 



L/2 



dxip^-\x) 



dx 

(37) 

d 2 ^ n +/ -\x) 



dx 2 



(38) 



It may be seen that (A 2 )qm 



(X 2 



CL 



' 1 ' O 

large n limit, in which case the uncertainty product 



i in the 



lim (AX) (aP) 

n->oo V / QM V / 



QM V /QM 

From the commutator relation, 

2 L 
L nirn 



(39) 



[X,P] 



2i 
n7r' 



it is clear that 



AX 



QM 



AP) 



> 



J QM n 2 7T 2 



(40) 



(41) 



and in the large n limit one recovers the expected result 
AX) (AP) > 0. However, it may be noted 

J QM \ J QM 

that the stationary state uncertainty product (|39[) does 
not vanish in the limit n — > oo - but it approaches the 
value ^ - which coincides exactly with that associated 
with the classical ensemble (see (|3"Tj) L 



C. Bouncing ball 

We now consider the example of a particle bouncing 
vertically up and down in a uniform gravitational field, 
which is described by the confining potential, 



V(z) = 



oo for z < 
mgz for z > 0. 



(42) 



A classical particle of mass m, energy E is under the 
influence of a constant force F — = —mg, z > 

and it bounces back and forth between < z < A, where 
A = E /mg denotes the maximum height reached. 

An ensemble of bouncing balls of energy E is character- 
ized by the classical position probability distribution Q 



-r^ for < z < A 

Pci J (z)=l 2A V^Ta) - - (43) 

I otherwise. 

which is readily identified by substituting (j4"2"]l in ([3])). 

Employing dimensionless position and momentum 
variables 



Z 



A' 



P 



V2mE y/2m 2 g A 



(44) 



(so that < Z < 1 and -1 < P < 1 for the 
bouncing particles) we obtain the classical moments 
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(Z)cl, (Z 2 )cl, (P)cl, (P 2 )cl as, 
(Z)cl = ~J f dzP C h(z)z 
^ dz ■ 



2 A 2 



,o y/1 - (z/A) 3' 
(Z 2 )cl = ^2 / dzP Ch (z)z 2 



i r A . 



(45) 



(46) 



2 v / 2m 2 g A Jo 



+ y / '2m(E — mgz) = 0, 



dz Pcl(z) (~y / 2m(E — mgz) 

(47) 



1 



>CL = 



2m 2 g A 



2A 



dz Pcl(z) 2m {E — mgz) 



1 



= — / dz y/1 ( Z /A) = - 



(48) 



Thus, the variances of Z and P are given by (AZ)^ L = 
and (AP)c L = |, leading to 



(AZ) 2 L (AP) 2 L 



135 



(49) 



Stationary state solutions of a quantum bouncer f]~Qj ] 
(confining potential of which is given by (T4"2"10 arc 
obtained by solving the time- independent Schrodinger 
equation 



ft 2 d 2 i, n {z) 



+ mgzip n (z) = E n ip n (z), (50) 



2m dz 2 
with the boundary condition 

^n(O) = 0. (51) 

In terms of the characteristic gravitational length 10] 



.2 \ 1/3 



1 2 m 2 g 



it is convenient to define dimensionless quantities 



EL 



i _ * ^/ 

mglg ' Zg " 



(52) 



(53) 



so that the Schordinger equation ((50)) takes the standard 
form 



d2Mz ' ] -z'ti (z') 
dz' 2 - z ^Az) 



(54) 



which is the Airy differential equation. The solutions of 
([54]) are two linearly independent sets of Airy functions 
Ai(z'), Bi(z'); however, the function Bi diverges as its 



argument increases, and so it is not a physically admissi- 
ble solution. The stationary state solutions of a quantum 
bouncer are thus given by, 



i> n {z') = N n Ai(z'), z'>-E' n 



(55) 



n 


E' n 


1 


2.3381 


2 


4.0879 


3 


5.5205 


4 


6.7867 


5 


7.9441 



TABLE I: The first few scaled energy eigenvalues E' n of quan- 
tum bouncing ball. 



where N n denotes the normalization constant. From the 
boundary condition (jSTj) . one obtains Ai(— E' n ) = 0, n — 
1,2,... leading to the identification that the (scaled) en- 
ergy eigenvalues E' n are the nth zeros of the Airy func- 
tion. The first few energy eigenvalues E' n (first few zeros 
of the Airy function) of the quantum bouncing ball are 
given in Table 1. 

Identifying the classical turning point A n associated 
with the energy eigenvalues E n of the quantum bouncer 
to be 



Ar 



= h E 'n 

mg 



(56) 



we define appropriately scaled position and momentum 
operators (which are quantum counterparts of Z, P de- 
fined in (011)) as, 



z = T 



IgK 



(57) 



P 



P 



\/2mE„ 



l gP 



Further, substituting (|5"2"]) . (|5"5)) in ([ST)) , we may express 
the configuration representation of the operators Z, P in 
terms of z' , E' as follows: 



P -> 



d 



rw dz 1 



(58) 
(59) 



The expectation values (Z) QM , (Z 2 ) QM , (P)qm, (P 2 )qm 
are then evaluated (numerical computation using Math- 
cmatica package) in the eigen states ((55)) of the quantum 
bouncing ball: 
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{Z)qm 



IQM 



EL 



N 2 f°° 

/ dz'A l 2 (z')(z' + E' n 

J —EL 



E' 2 



dz' \^ n {z')\ 2 {z 1 + E' n f 



N 2 f°° 8 

^ J dz'Ai>(z')(z' + E' n f = ± 



(60) 



(61) 



(P)qm 



(P 2 )qm 



VK J -EL. 

-iN 2 r°° 

J-E: 



dz lj)*(z ) 



, dlp n (z') 



dz' 



i 



dAi(z ) 
dz' Ai(z') — t 1 -^- 
dz' 



= o 



(62) 



(63) 



The quantum expectation values (|60l) - (j6"3")l in the eigen 
states match identically with those (see ([4"5"j) -(ji5 |) ) asso- 
ciated with a classical ensemble of bouncing balls. This 
is a novel identification, bringing forth the deep rooted 
unifying features in classical and quantum realms. 

From Eqs. (1501) - (| 6 3 1) , we obtain the variances of Z 
and P in the stationary states to be (AZ)qm = 4/45, 
(AP)qm = 1/3 and hence, the uncertainty product 

(AZ) 2 QM (AP) 2 M = 4-, (64) 



135' 



which matches exactly with that of the classical ensemble 
of bouncing balls (see (03)). It may be noted that, the 
commutation relation 



[Z,P] 




(E' n f/ 2 



(65) 



would lead to the uncertainty relation 

( A ^)qm ( AP )qm > TiWJt- In the lar S e n limit 
— > (as the energy eigenvalues obey the scaling rela- 
tion Q E' n oc n 2 / 3 with n), thus resulting in the classical 
limit on the variance product (AZ)q M (AF)q M > 0. 
The exact result on the uncertainties (see ([M])'). which 
holds for all the eigenstates of the quantum bouncing 
ball, approaches the value j|g for all values of n and it 
agrees perfectly with that associated with the classical 
ensemble. 

IV. SUMMARY 



h — > and also, to that based on Eherenfest's theorem. 
However, both these descriptions on classical regime are 
found to be unsatisfactory Quantum predictions 

being purely statistical in nature, it is pertinent to com- 
pare the classical limit of the quantum system with the 
corresponding classical ensemble - not with a single parti- 
cle 0, IH . And a suitable operational criterion is therefore 
to compare the averages and probability distributions in 
both realms. In this connection, we have explored paral- 
lels between variances of dimensionless position and mo- 
mentum variables - and hence the resulting uncertainty 
relations - in the stationary eigen states of one dimen- 
sional bound quantum systems and their classical coun- 
terparts. In the specific examples of harmonic oscillator, 
infinite well and the bouncing ball, we have shown that 
the first and second moments of the scaled canonical ob- 
servables exactly agree with each other - resulting in iden- 
tical uncertainty relations in both quantum and classical 
regimes. This identification reflects a deep underlying 
connectivity between the two formalisms - despite their 
differences in the mathematical and conceptual basis. 
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